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Recent more precise experiments have revealed both a slow and a fast quantum oscillation in the c-axis 
resistivity of nearly optimal to overdoped electron-doped high temperature superconductor Nd2- x Ce x Cu04. 
Here we study this problem from the perspective of Fermi surface reconstruction using an exact transfer matrix 
method and the Pichard-Landauer formula. In this method, neither quasiclassical approximations for magnetic 
breakdown, nor ad hoc broadening of Landau levels, are necessary to study the high field quantum oscillations. 
The underlying Hamiltonian is a mean field Hamiltonian that incorporates a two-fold commensurate Fermi 
surface reconsruction. While the specific mean field considered is the d-density wave, similar results can also 
be obtained by a model of a spin density wave, as was explicitly demonstrated earlier. The results are consistent 
with an interplay of magnetic breakdown across small gaps in the reconstructed Fermi surface and Shubnikov-de 
Haas oscillations. 



I. INTRODUCTION 

Quantum oscillations were first discovered 1 in the Hall 
coefficient of hole-doped high temperature superconductor 
YBa2Cu306+5 (YBCO) at high magnetic fields between 
35 — 62T in the underdoped regime close to 10%. Since then 
a number of measurements, in even higher fields and with 
greater precision using a variety of measurement techniques 
have confirmed the basic features of this experiment. How- 
ever, the precise mechanism responsible for oscillations has 
become controversial. 2 Fermi surface reconstruction due to a 
density wave order that could arise if superconductivity is "ef- 
fectively destroyed" by high magnetic fields have been focus 
of some attention!^ 

In contrast, similar quantum oscillation measurements in 
the doping range 15 - 17% in Nd 2 _ x Ce x Cu0 4 (NCCOj^ 
seem easier to interpret, as the magnetic field range 30 — 65T 
is far above the upper critical field, which is less than 10T. 
This clearly places the material in the "normal" state, a source 
of contention in measurements in YBCO; in NCCO the crystal 
structure consists of a single CuO plane per unit cell, and, in 
contrast to YBCO, there are no complicating chains, bilayers, 
ortho-II potential, stripes, etcP Thus, it would appear to be 
ideal for gleaning the mechanism of quantum oscillations. On 
the other hand, disorder in NCCO is significant. It is believed 
that well-ordered chain materials of YBCO contain much less 
disorder by comparison. 

In a previous publication!^ we mentioned in passing that 
it is not possible to understand the full picture in NCCO 
without magnetic breakdown effects, since the gaps are ex- 
pected to be very small in the relevant regime of the parame- 
ter space. However, in that preliminary work the breakdown 
phenomenon was not addressed; instead we focused our at- 
tention to the effect of disorder. Since then recent measure- 
ments^ have indeed revealed magnetic breakdown in the range 
16 — 17% doping, almost to the edge of the superconduct- 
ing dome. Here we consider the same transfer matrix method 
used previously, 6 but include third neighbor hopping of elec- 
trons on the square planar lattice, without which many ex- 
perimental aspects cannot be faithfully reproduced, including 
quantitative estimates of the oscillation frequencies and break- 



down effects. The third neighbor hopping makes the numer- 
ical transfer matrix calculation more intensive because of the 
enlarged size of the matrix, but we were able to overcome the 
technical challenge. In this paper we also analyze the c-axis 
resistivity and the absence of the electron pockets in the ex- 
perimental regime. 

n. HAMILTONIAN 

The mean field Hamiltonian for <i-density wav tP(DDW) in 
real space, in terms of the site-based fermion annihilation and 
creation operators c\ and q, is 

Hddw = X/ iC i Ci + eiai jc i c J + h - c -> (!) 

where the nearest neighbor hopping matrix elements include 
DDW gap Wq and are 



where (n, m) are a pair of integers labeling a site: i = nx + 
my; the lattice constant a will be set to unity unless otherwise 
specified . In this paper we also include both next nearest hop- 
ping matrix element, t', and third nearest neighbor hopping 
matrix element t". A constant perpendicular magnetic field 
B is included via the Peierls phase factor ai j = ^ J? A • dl, 
where A ~ (0, —Bx, 0) is the vector potential in the Landau 
gauge. The band parameters are chosen to be t — 0.38eV, 
t' = 0.32i, and t" = OM'W The chemical potential fi is 
adjusted to achieve the required doping level and is given in 
Table [I] so is the DDW gap Wq. We assume that the on-site 
energy is (5-correlated white noise defined by the disorder av- 
erage e[ = and e[ej = V^ 2 <5i j . Disorder levels for each of 
the cases studied are also given there in Table [I] We have seen 
previously that longer ranged correlated disorder lead to very 
similar results.^ 
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FIG. 1. (Color online) A plot showing the breakdown junctions and 
electron trajectories across them in the extended Brillouin zone. The 
figure corresponds NCCO with 17% doping and small DDW gap. 
Note that the reflection at the junctions involve a large change in mo- 
mentum. The electron trajectories that lead to magnetic breakdown 
of small hole pockets are shown. 



The Fermi surface areas (See Fig. [TJ of the small hole 



pocket in the absence of disorder correspond to oscillation 
frequencies 330T at 15% doping, 317T at 16% doping and 
291T at 17% doping. These frequencies seem to be insensi- 
tive to Wq within the range given in Table [I] 



III. THE TRANSFER MATRIX METHOD 

Transfer matrix to compute the oscillations of the conduc- 
tance is a powerful method. It requires neither quasiclassical 
approximation to investigate magnetic breakdown nor does it 
require ad hoc broadening of the Landau level to incorpo- 
rate the effect of disorder. Various models of disorder, both 
long and short-ranged, can be studied ab initio. The mean 
field Hamiltonian, being a quadratic non-interacting Hamil- 
tonian, leads to a Schrodinger equation for the site ampli- 
tudes, which is then recast in the form of a transfer matrix; 
the full derivation is given in the Appendix. The conductance 
is then calculated by a formula that is well known in the area 
of mesoscopic physics, the Pichard-Landauer formula.!^ This 
yields Shubnikov-de Haas oscillations of the a&-plane resis- 
tivity, p a t,. We show later how this can be related to the c-axis 
resistivity p c measured in experiments. 

Consider a quasi- ID system, N M, with a periodic 
boundary condition along y-direction. Here Na is the length 
in the ^-direction and Ma is the length in the y-direction, a 
being the lattice spacing. Let <J„ = (^„,i, V>n,2, ■ • ■ , ^n.Ai) T , 
n = 1 , . . . N, be the amplitudes on the slice n for an eigen- 
state with a given energy. Then the amplitudes on four suc- 
cessive slices must satisfy the relation 
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where U nj A nj B n , C nj D n are M x M matrices. The non- 
zero matrix elements of matrix A n are 



(-^n)m, ?7i+l 
(^n)m,m— 1 — 



i 



(-1) 



m+n 



—i'eK-"-!)^ 
_£' e i(n+i)<A 



(4) 

(5) 
(6) 



where <f> = Ba 2 e/hc is a constant. The elements of the matrix 



B„ are 
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t"l, where 1 is the M x M 



Here C n = A\ x and D n 
identity matrix. 

The AM Lyapunov exponents, 7.;, of limjv_ J . 00 (T/v7^), 

where Tn — Y[ J j=i Tj, are defined by the corresponding 
eigenvalues = e 7i . All the Lyapunov exponents 7! > 
72 > ... > 74M, are computed by a method described in 
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TABLE I. Parameters Wo (DDW gap), Vo (on-site disorder poten- 
tial), and fi (chemical potential). 
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Ref. [121 However, the matrix is not symplectic. Therefore 
all AM eigenvalues are computed. Remarkably, except for a 
small set, consisting of large eigenvalues, the rest of the eigen- 
values do come in pairs (A, 1/A), as for the symplectic case, 
within our numerical accuracy. We have no analytical proof 
of this curious fact. Clearly, large eigenvalues contribute in- 
significantly to the Pichard-Landauei 1 1 formula for the con- 
ductance, a a b(B): 



2M 



1 (TnU) 



(12) 



We have chosen M to be 32, smaller than our previous work. 6 
The reason for this is that the matrix size including the third 
neighbor hopping is larger AM x AM instead of 2M x 2M. 
We chose N to be of the order of 10 6 , as before. This easily 
led to an accuracy better than 5% for the smallest Lyapunov 
exponent, 7^ in all cases. 



IV. MAGNETIC BREAKDOWN AND QUANTUM 
OSCILLATIONS 



We compute the conductance as a function of the magnetic 
field and then Fourier transform the numerical data. This pro- 
cedure of course depends on the number of data points sam- 
pled within a fixed range of the magnetic field, typically be- 
tween 45 — 60T. As the number of sampling points increases, 
the peaks become narrower but greater in intensity, conserv- 
ing the area under the peak. But the location of each peak 
and the relative ratio of the intensities remain the same. In 
order to compare the Fourier transformed results, we keep the 
sampling points fixed in all cases to be 1200. 

In Figs.|2]through|4]the results for 17% doping for a 5meV 
gap and varying degrees of disorder are shown. Both the 
slow oscillation at a frequency 290T corresponding to the 
small hole pocket and 11, 700T corresponding to the large 
hole pocket, as schematically sketched in Fig. [T] in the ex- 
tended Brillouin zone, can be seen. Note that partitioning 
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FIG. 2. (Color online) Fourier transform of the conductance oscil- 
lations with a smooth background term subtracted. The parameters 
correspond to 17% doping with a DDW gap of 5meV and disorder 
Vo = 0.2*. The horizontal axis is in units of Tesla and the vertical 
axis is in arbitrary units. 
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FIG. 3. (Color online) The same as in Fig. [2] corresponding to 17% 
doping but with a DDW gap of 5meV and disorder Vo = 0.4*. 



of the spectral weight between the peaks changes as the de- 
gree of disorder is increased. If we change the value of the 
gap to lOmeV, shown in Figs. [5] through [7] the overall pic- 
ture remains the same, although the slower frequency peak 
is a bit more dominant, as the magnetic breakdown is a lit- 
tle less probable. For 16% doping similar calculation with 
gaps of lbmeV and 30meV also show some evidence of mag- 
netic breakdown depending on the disorder level, particularly 
seen in IhmeV data in Figs. [8]though 10 On the other hand, 
the evidence of magnetic breakdown is much weaker in the 
30mel^ data shown in Figs. [TT| through 13 



It is important to note that in none of these calculations one 
finds any evidence of the electron pocket centered at (tt, 0) 
and its symmetry counterparts, which should roughly corre- 
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FIG. 4. (Color online) The same as in Fig. [2] corresponding to 17% 
doping but with a DDW gap of 5meV and disorder Vo = 0.6*. 
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FIG. 5. (Color online) The same as in Fig. [2] corresponding to 17% 
doping but with a DDW gap of WmeV and disorder Vo = 0.2t. 
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FIG. 8. (Color online) The same as in Fig. [^corresponding to 16% 
doping but with a DDW gap of 15meV and disorder Vo = 0.2t. 
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FIG. 6. (Color online) The same as in Fig. [2] corresponding to 17% 
doping but with a DDW gap of WmeV and disorder Vo ~ OAt. 
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FIG. 9. (Color online) The same as in Fig. [2] corresponding to 16% 
doping but with a DDW gap of 15meV and disorder Vo = OAt. 



spond to a frequency of 2700T. This is in part due to the fact 
that the effect of disorder is stronger on the electron pocket 6 
and in part due to the fact that at the breakdown junctions 
transmission coefficient is larger than the reflection coefficient 
because it entails a large (tt/2) change in the direction of the 
momentum; see Fig. [T] 



V. OSCILLATIONS IN c-AXIS RESISTIVITY 

The Pichard-Landauer formula was calculated for conduc- 
tance oscillations in the afe-plane, while the actual measure- 
ments in NCCO are carried for the c-axis resistivity. It is 
therefore necessary to relate the two to compare with experi- 
ments. A simple description for a strongly layered material 
can be obtained by modifying an argument of Kumar and 
Jayannavar. 13 An applied electric field, E, along the direction 



perpendicular to the planes will result in a chemical potential 
difference 



A/i = edE, 



(13) 



where d is the distance between the two planes of an unit cell. 
The corresponding current, j c , is (e p is the Fermi energy) 



j c = e [A/j,g 2D (e F ,H)] 7, 



(14) 



since A/j,g 2 D(£F, H) is the number of unoccupied states to 
which an electron can scatter, while 7 is the scattering rate 
between the planes of a unit cell. Here, we have included 
a possible oscillatory dependence of the the two-dimensional 
density of states, g 2 D(£F, H), that gives rise to Shubnikov-de 
Haas oscillations in the afr-plane. Thus, 
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FIG. 7. (Color online) The same as in Fig. [2] corresponding to 17% 
doping but with a DDW gap of WmeV and disorder Vo = 0.6t. 



FIG. 10. (Color online) The same as in Fig. [^corresponding to 16% 
doping but with a DDW gap of 15meV and disorder Vo = 0.6t. 
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FIG. 11. (Color online) The same as in Fig. [^corresponding to 16% 
doping but with a DDW gap of 30meV, and disorder Vo = 0.2t. 
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FIG. 12. (Color online) The same as in Fig. [^corresponding to 16% 
doping but with a DDW gap of 30meV and disorder Vo ~ OAt. 



There is an implicit assumption: an electron from a given 
plane makes a transition to a continuum of available states 
with a finite density at the Fermi surface in the next plane. We 
are not interested in the Rabi oscillations between two discrete 
states, a process that cannot lead to resistivity. 

The measured a6-plane resistivity is of the order 10/iil-cm 
as compared f2-cm for the c-axis resistivity even at optimum 
doping, 5 which allows us to make an adiabatic approximation. 
Because an electron spends much of its time in the plane, mak- 
ing only infrequent hops between the planes, we can adiabati- 
cally decouple these two processes. The slower motion along 
the c-axis can be formulated in terms of a 2 x 2 matrix for 
each parallel wave vector kn after integrating out the planar 
modes. For simplicity, we are assuming that the c-axis warp- 
ing is negligible, so there are only two available states of the 
electron corresponding to its locations in the two planes. The 
excitations in a plane close to the Fermi surface, kn s» k F n, 
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FIG. 13. (Color online) The same as in Fig. [^corresponding to 16% 
doping but with a DDW gap of 30meV and disorder Vo = 0.6t. 



can be approximated by a bosonic heat bath of particle-hole 
excitations. In this language, the problem maps on to a two- 
state Hamiltonian 



H — —t c a x 



E 



fujjb]bj 



Mb 



(16) 



where cr's are the standard Pauli matrices and t c is the hopping 
matrix element between the nearest neighbor planes.. Given 
the simplification, the sum over k|| is superfluous, and the 
problem then maps on to a much studied model of a two-level 
system coupled to an Ohmic heat bath. ^ The Ohmic nature 
follows from the fermionic nature of the bath. 15 The effect of 
the bath on the transition between the planest is summarized 
by a spectral function, 



J{ U ) = ~'£f*5(u,-u, j ). 



For a fermionic bath, we can choose 
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where uj c is a high frequency cutoff, which is of the order of 
uj c = 2 /T a b, where r ao is of the order of the planar relax- 
ation time. For a Fermi bath, the parameter a is necessarily 
restricted to the range < a < 1.^ Moreover, for coherent 
oscillations we must have a < 1/2.™ However, we shall leave 
a as an adjustable parameter, presumably less than or equal to 
1 /2 to be consistent with our initial assumptions. While a sim- 
ilar treatment is possible for a non-Fermi liquid, 16 the present 
discussion is entirely within the Fermi liquid theory. 

The quantity 7 is the interplanar tunneling tunneling rate 
renormalized by the particle-hole excitations close to the pla- 
nar Fermi surface and can be easily seen to be^ 



2t c 

The c-axis resistivity is then 
h 1 



2tc 

hid. 



e 2 dg 2 D(£F, H)hhj c 



Hu> c 
2t~r. 



(19) 



(20) 



This equation can be further simplified by expressing it as a 
ratio of p c j p ao , but this is unnecessary. Two important qual- 
itative points are: p c is far greater than p a t, and the root of 
the quantum oscillations of p c is quantum oscillations of the 
planar density of states. 



VI. CONCLUSIONS 

We have shown that a qualitatively consistent physical pic- 
ture for quantum oscillations can be provided with a simple set 
of assumptions involving reconstruction of the Fermi surface 
due to density wave order. Although the specific order con- 
sidered here was the DDW, we have shown previously that at 
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the mean field level a very similar picture can be provided by 
a two-fold commensurate spin density wave (SDW)!^ Thus, 
it appeared unnecessary to repeat the same calculations using 
the SDW order. 

In YBCO, studies involving tilted field seems to rule out 
triplet order parameter, hence SDW. 17 Moreover, from NMR 
measurements at high fields, there appears to be no evidence 
of a static spin density wave order in YBCO .^Similarly there 
is no evidence of SDW order in fields as high as 23. 2T in 
YBa 2 Cu40^51 while quantum oscillations are clearly ob- 
served in this materialPQf Also no such evidence of SDW 
is found up to 44T in Bi2Sr2_ x La x Cu06+,5p^ At present, 
results from high field NMR in NCCO does not exist, but 
measurements are in progress.^ It is unlikely that such static 
SDW order will be revealed in these measurements. This con- 
jecture is based on the zero field neutron scattering measure- 
ments which indicate very small spin-spin correlation length 
in the relevant doping regimeP^A long range SDW order can- 
not appear merely by applying high magnetic fields, which is 
energetically a weak perturbation even for 45T fieldP^l 

As to singlet order, most likely relevant to the observation 
of quantum oscillations] 2 ^ charge density wave is a possibil- 
ity, which has recently found some support in the high field 
NMR measurements in YBCO. 18 But since the mechanism 
involves oxygen chains, it is unlikely that the corresponding 
NMR measurements in NCCO will find such a charge order. 
As to singlet DDW, there are two neutron scattering mea- 
surements that seem to provide evidence for itPSl However, 
these measurements have not been confirmed by further inde- 
pendent experiments. However, DDW order should be con- 
siderably hidden in NMR involving nuclei at high symmetry 
points, because the orbital currents should cancel. 

A mysterious feature of quantum oscillations in YBCO is 



the fact that only one type of Fermi pockets are observed. If 
two-fold commensurate density wave is the mechanism, this 
will violate the Luttinger sum rulePEzl We have previously 
provided an explanation for this phenomenon in terms of dis- 
order arising from both defects and vortex scattering in the 
vortex liquid phased however, the arguments are not unas- 
sailable. In contrast, for NCCO, the experimental results are 
quite consistent with the simple theory discussed above. We 
have not addressed AMRO in NCCO, as the data seem to 
be somewhat anomalous, 7 although within the Fermi liquid 
framework discussed here, it should be possible to address 
this effect in the future. 

The basic question as to why Fermi liquid concepts should 
apply remains an important unsolved mystery. 28 It is possible 
that if the state revealed by applying a high magnetic field has 
a broken symmetry with an order parameter (hence a gap), the 
low energy excitations will be quasiparticle-like, not a spec- 
tra with a branch cut, as in variously proposed strange metal 
phases. In this respect, the notion of a hidden Fermi liquid 
may be relevant.^ 
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Appendix: The derivation of the transfer matrix 

The DDW Hamiltonian in real space is 



<iJ>" 



iW , 



(A.l) 



Here, e la ' j is the Peierls phase due to the magnetic field. The 
summation notations are as follows: (i, j)', and (i, j)" 

imply sum over nearest-neighbor, next-nearest-neighbor, and 
the third-nearest-neighobor sites, respectively. For example, 
with the lattice constant set to unity, is satisfied when 
i = j ± x or i = j ± y. Likewise, (i, j)' requires i = j + x ± y 
ori = j — x ± y and requires i = j ± 2x or i = 



j ± 2y. Here Wq is the DDW gap and i = (n, to). Consider 
an eigen state |^) with an energy eigenvalue E: H\^f) = 
E\^>), where 1^) = ^/>(i)|i); the amplitude at a site is 
tp(i). Then the Schrodinger equation can be written in terms 
of the amplitudes ijj n (m) of the n-th slice for all values of 
to = 1,2, . ..M: 
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(A.2) 



With periodic boundary condition along the y-axis i.e. 
tp n (M + 1) = ip n (l), the Schrodinger equation can be ex- 
pressed as a matrix equation: 

= -U n 1p n+ 2 + A n ip n+1 + B n ip n + C n 1p n -1 + D n 1p n _2, 

(A3) 

where U n , A n , B n , C n , and D n are M x M matri- 
ces defined in equations following Eq. [3] Now we can 



solve the Schrodinger equation for tpn+2 to obtain V'n+2 = 
t/,7 1 (A n *fj n+1 + B n ip n + C„Vn-i + AiV'n-a)- Then the am- 
plitudes at a set of four successive slices, ipn~i through ip n +2, 
can be written in terms of the amplitudes of a previous set of 
four successive slices, ip n -2 through ip n +i- Thus, the transfer 
matrix in the main text follows. 
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